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■ Abstract 

Q^' We show how radiative QCD corrections calculated in terms of quarks can be 

O . incorporated at the hadron level in inclusive semileptonic i3-meson decays. The 

bound state effects are described by a momentum distribution function of the b 
quark. The summation over the final states and the averaging over the momentum 
^ . distribution of the decaying quark render the radiative corrections finite. With this 

^ I coherent formalism we investigate the shape of the electron spectra for b ^ u and 

b ^ c decays as a function of the parameters of the theory. The resultant b ^ c 
electron energy spectrum is in agreement with the experimental data. 
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1 Introduction 



Semileptonic i?-meson decays have been studied for some time now. The inclusive decays 

B ^ Xg + e + Ue, (1) 

with B representing 5^ or B^ and Xg any possible hadronic final state containing a charm 
quark (g = c) or an up quark {q = u), provide information on both couplings Vcb and Vub, 
as well as new information on the internal structure of the i?-meson. In this paper we 
study the electron energy spectra. For the theoretical description of the electron spectra 
in eq.(l), strong interactions in the underlying weak decays must be incorporated, since 
they are responsible for the confinement of quarks and gluons into hadrons. They will 
be included in two steps: as bound state effects and also in the form of gluons radiated 
during the decay. 

For inclusive i?-meson decays, it was recognized that extended regions of phase space 
involve large values of g^, the momentum transfer squared, originating from the large mass 
of the B meson. Consequently the commutator of the two currents describing the decay 
is dominated by distances close to the light-cone. In this case it is justified to replace the 
commutator of the weak currents by their light-cone singularity times a bilocal operator 
in 6-quark fields [1-3]. This replacement together with standard mathematical methods 
leads to general expressions of the decay spectra which involve a 6-quark distribution 
function, whose origin is non-perturbative. Thus the semileptonic decays in eq. (|I]) can be 
described, in direct analogy to deep inelastic scattering, in terms of a quark distribution 
function, which depends on a new scaling variable 

It has also been recognized that the heavy quark field can be studied in an effective field 
theory derived from the QCD Lagrangian. The heavy quark effective theory (HQET) sets 
a framework for keeping track of the heavy quark mass dependence and for parametrizing 
nonperturbative phenomena. The effective theory has been successfully applied to inclu- 
sive B decays [4-12]. Using the operator product expansion and the method [5-7] of the 
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HQET, it was possible to derive sum rules for the distribution function, which depend on 
the kinetic energy and the chromomagnetic energy of the 6-quark in the i?-meson. The 
numerical values for the sum rules are determined by static properties of B mesons and 
QCD sum rules |T^. The sum rules specify the mean value and the variance of the dis- 
tribution function. They imply that the distribution function /(^) peaks at large values 
of ^ ~ 0.93 and is very narrow. 

In addition radiative QCD corrections must be included. The QCD radiative correc- 
tions to the electron energy spectra were computed [14-18] at the quark level. In the 
quark decay the phase space ends at the electron energy Ee = ^ — , with nib and 
rriq the masses of the 6-quark and the final quark, respectively. The quark decay rates 
as well as the radiative corrections vanish for E^. > ^(1 — ^), whereas in reality the 



b 
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physical endpoint is Ee = ^{1 jf¥^) with M the mass of the 5-meson and Mx^i„ the 

minimum value of the invariant mass of the hadronic final state. In addition, the 0{as) 
radiative corrections at the quark level have logarithmic singularities at the endpoints. 
For inclusive decays, however, we integrate over the phase space of the final quark and 
average over the momentum distribution of the initial quark in order to incorporate the 
bound state effect. These steps sum over ensembles of states, render the radiative cor- 



rections finite ||19| , pO| and extend the phase space from the quark level to the hadron 
level. 

In this way we have at our disposal a coherent treatment of perturbative and non- 
perturbative QCD effects. This treatment has the advantage of accounting correctly for 
the phase-space effects and producing the electron energy spectra which are smooth ev- 
erywhere up to the physical endpoints. For instance, the spectrum is a smooth func- 
tion of the electron energy Eg in the endpoint region between the b ^ c endpoint 
E^ = = 2.31 GeV and the 6 ^ m endpoint E^ = f = 2.64 GeV, which 

is useful for extracting \Vub\ since only the b u transition is allowed in this region. 

For the sake of completeness we present in section 2 the general formalism for inclusive 
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semileptonic B meson decays. Then we show in section 3 that the hght-cone dominance 
leads to an expression for the decay rate in terms of a distribution function. In these two 
sections, we point out how the general formula in eq.(7) reduces to eq.(24). Properties 
of the distribution function are discussed in section 4 where they are quantified by the 
heavy quark effective theory. In section 5 we propose a method for including the QCD 
radiative corrections to bound states. With the above formalism we calculate and analyse, 
in section 6, the electron energy spectra. 



2 General Formalism 

The inclusive semileptonic decays (|I|) , in which the B meson of four-momentum P decays 
into an electron of four- momentum kg and an antineutrino of four- momentum k,^, are 
described by the decay amplitude 

M = V,,^u{k,)Y{l - l,Hku){n\j,mB) . (2) 

Here Vgb are the elements of the CKM matrix and j^{x) is the weak current, which in 
terms of quark fields is given by 

= Q{xh^^{l - l5)b{x) (3) 

and l-B) is the B-meson state normalized according to {B\B) = 2 Pq (271)^6^(0). The basic 
quantity for the decay is the second rank tensor 

- [27Trd-[P-q- 

i=l 



n 



\ (27r)32Ei_ 



{27:f5\P - q -Y.mB\m\n){n\3,mB) , (4) 

i=l 

where q stands for the four-momentum transferred from the decaying B meson to the 
lepton pair, q = k(. + k^. It is useful to express the hadronic tensor in terms of a current 
commutator 

W,. = -^l d'ye^'^-y{B\[j,{y),jl{0)]\B) (5) 

because the commutator is more convenient for theoretical considerations. The hadronic 
tensor can be decomposed in terms of scalars Wa{q'^, q ■ P), follows : 

W,. = -g,.W, + - .W^^3 + ^W^4 + ^^^^^^P^W, . (6) 



The tensor {P^qi, — q^^Pu) does not appear because of the time reversal invariance. We 
can express the differential decay rates in terms of the five hadronic structure functions 
Wa, a = 1, . . . , 5. The decay rate of the process (|ID in the rest frame of the B meson is 



Gp\Vqh\'^ 



2 2 

W,q^ + W2{2E,qo - 2EI - |-) + W^^{q^ 



2K 



(7) 



The structure functions and do not appear above because their contribution is 
proportional to the square of the electron mass and we ignore the lepton masses. In this 
general formalism the unknown hadronic structure resides in the functions Wa- 



3 Light-Cone Dominance 



It is well known that integrals like the one in eq.(^ are dominated by distances where 

0<y'<\. (8) 
q^ 

For inclusive semileptonic B-meson decays (|l|), is timelike and varies in the physical 
range 

0<q'<iM-Mx^J'. (9) 

For extended regions of phase space the momentum transfer squared satisfies > g^g^ 
with g^gj ~ 1 GeV^. In these regions we expect the dominant contribution to the integral 
in eq.(^ to come from distances of the current commutator close to the light-cone. The 
commutator in this region is in fact singular leading to the dominant contribution 



{B\ [j,{y),jl{0)\ \B) = 2{S,^,p - ie,^,p) [9"A,(y)] (5|6(0)7^(1 - l,)h{y)\B) , (10) 
where S^^^p = g^aQvp + g^pQua - Qi^uQap and Ag(?/) is the Pauli- Jordan function for a free 



g-quark of mass m^. The factor in the square bracket in eq.([TO|) with the derivative of 
the Pauli- Jordan function has a singularity on the light-cone. The last factor with the 
reduced matrix element contains the long-distance contribution. The product of those 
two factors is Lorentz covariant and can be calculated in any Lorentz frame of reference. 
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The reduced matrix element has a simple Lorentz structure. It is in general a function 
of two scalars y"^ and y ■ P and can be expanded in powers of y"^: 

oo 

(5|6(0)7^(i - i,)b{y)\B) = 47rP^ T^iv'TPniy ■ P) ■ (n) 

n=0 

We shall keep the first term of the series because the higher order terms are suppressed 
by powers of q~^. This approximation is justified provided the coefficients are not very 
large. We can estimate the coefficients in quark models or the heavy quark effective theory 
which indicate that they satisfy 

{y'TF^iy . P) ^ e-'^^^-y{Klcn/q'T , (12) 

where v is the velocity of the initial B meson, defined by f = P/M. This behavior 
motivates the truncation of the series by keeping the first term with n = 0. 

The Fourier transform of -fo(l/ ■ P) defines the quark distribution function 

/(O = -r^, I d{y.P)e^^y^{Bmm-l^)Ky)\B)\y^=,. (13) 



47rM2 

We can use the inverse Fourier transform 



F,{y ■P) = ^ I die-^^y^'fiO (14) 



and substitute Fq in eg s. ([Til), ( [T0|) and (||), then carry out the y-integration in eq.(|) and 
arrive at 



4(5^«.;3 - ie,o,,p) J d^fiO^m - qo)S [{^P - qf - m^] (^P - q^P^ . (15) 



The components of the tensor W^^ are expressed in terms of the distribution function. 
We have shown that the dominance of the light-cone makes possible the expression of the 
decay rate in terms of a quark distribution function defined in eq.(13). 



A special consequence of the decay kinematics is the occurrence of two roots in the 
argument of the 5-function in eq.(|T5|), namely 



q- P± J{q- Py - M^{q^ - ml) 
= M'2 • (16) 
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We shall elaborate on this property below. The light-cone dominance ascribes the five 
hadronic structure functions to a single light-cone distribution function. The explicit 
relations are the following 



W, = 2[/(e+)+/(e-)], (17) 

W2 = -^[e+/(e+)-e-/(e-)], (18) 

Ws = --^[/(e+)-/(e-)], (19) 

W4 = 0, (20) 

W5 = W3. (21) 



The structure functions are evaluated in two variables The second root, is 
a straightforward consequence of the analysis and corresponds to the creation of quark- 
antiquark pairs through the Z-diagram because the energy of the final quark is negative. 
The kinematic ranges for C,± are 

?r < e+<l, (22) 
^ < i-<l ^. 23 

In the light-cone and away from the resonance region /(^-) is relatively small. For b ^ c 
decays ^_ ~ 0.75 where /(^-) is negligibly small. Scaling of the structure functions with 
the scaling variable ^+ holds when /(^-) is negligible P]. 

The expression of the structure functions Wa in terms of a single distribution func- 
tion, which depends on two values ^± of the scaling variable, is a large simplification. 
Substituting the structure functions in eq.(0) we arrive at 



'° ^\ {m^m^E^M - q^) - (e+ - e-)} . (24) 



dEedq^dqo Atx^M /q2 + ml 



The remaining unknown function is the reduced matrix element on the light-cone 
whose Fourier transform appears as the 6-quark distribution function. 
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4 Properties of the Distribution Function 



The distribution function obeys positivity and is zero for ^ < ot ^ > 1 Three sum 
rules for the 6-quark distribution function are known. The first one expresses the b quark 
number conservation |3ll 



(25) 



Performing the operator product expansion to reduce the bilocal operator to local ones 
and following [5-7] in order to expand the matrix element of the local operator in the 
HQET, two more sum rules were derived. They determine up to order {^Qco/^bY ^^e 
mean value /x and the variance of the distribution function, which characterize the 
position of the maximum and its width, respectively: 



a 



Ef 



where 



2M 



B 



Ami 



2M 



B 



" 2ml 



B 



B 



(26) 
(27) 

(28) 
(29) 



with Eb = Gb + Kb- The first matrix element Gb parametrizes the chromomagnetic energy 
arising from the b quark spin and is determined by the mass splitting between B* and B 
mesons [5-7]. For the observed difference Mb* — Mb = 0.046 GeV 



mbGb = --{Mb* - Mb) = -0.034 GeV. 



(30) 



The second matrix element Kb parametrizes the kinetic energy of the b quark in the B 
meson. It is determined with the help of QCD sum rules and carries a larger error, leading 



to the result 21 



2mlKb = 0.5 ± 0.2 GeV^ 



(31) 
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Taking = 4.9 ± 0.2 GeV, the mean value and the variance of the distribution function 
are estimated to be 

/i = 0.93 ±0.04, (32) 

= 0.006 ±0.002, (33) 

indicating that the distribution function is sharply peaked around its mean value, which 
is close to one. These results are consistent with the original expectations that the distri- 
bution function of a heavy quark peaks at a large value of its argument. 



5 QCD Radiative Corrections 

Special attention must be paid to the radiative corrections from the emission of gluons 
and the associated virtual diagrams. QCD radiative corrections to the electron energy 
spectra were studied in several articles [14-18], where they were calculated at the quark 
level. As already mentioned in the introduction, in the application of radiative corrections 
at the hadron level we encounter two problems: the first is to change the quark phase 
space to the physical one and the second is the treatment of the logarithmic singularities 
to order a^, which appear at the quark-level endpoints = ^(1 — ^). 



These problems may be solved by taking into account the bound state effect. In the 
decay of a 5-meson, the perturbative QCD correction will be modified by the bound state 
effects, since QCD confinement implies that free quarks are not asymptotic states of the 
theory. The bound state effect is described by the 6-quark distribution function given in 



GQ-dUD) which is the probability of finding a b quark with momentum inside the B- 
meson. The substitution of the b quark momentum p^hj introduces the hadronic phase 
space. Furthermore, the radiative corrections obtained perturbatively must be convoluted 
with the distribution function. The final contribution for the radiative corrections is given 
by 

L^JUioj^^] , (34) 
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where the quark-level 0{as) perturbative QCD correction dT'^^^/dEe was computed an- 
alytically in p!7| , . In this way the endpoints of the perturbative spectra are extented 
from the quark level to the hadron level and the logarithmic singularities are eliminated. 
As we shall see, the interplay between perturbative and nonperturbative QCD effects 
is important, especially for the shape of the h ^ u electron energy spectrum near the 
endpoint. 

We could give at this point formulas for the radiative corrections at the quark level. 



but since they are available in two articles ||T^, [T^ we refer to them. The interested reader 
may consult these articles and use their formulas for dT^^^/dEe which occurs in eq.(34). 
An important property of eq.(34) is that the integral over ^ eliminates the logarithmic 
singularities. 

In order to calculate the decay spectra we need a distribution function f{^) consistent 
with the properties of section 4. We propose the Ansatz 

m = Nj^^^^^e{oe{i-0, (35) 

where N is the normalization constant and a and b two parameters. In case a = 
and b = nih/M, this distribution function reduces to a delta function, S{C, — nih/M), 
and thus reproduces the free-quark decay model. In addition, the constraints ( ^2l) and 
(B3), stemming from the HQET, limit the two constants a and b. Other forms of the 



distribution function have been proposed in p2|, p3|. 



We use eqs.(^4D, (|3^) and (|35D to compute the electron energy spectra for 6 — > c and 
b ^ u decays. The calculation is done using (|^) by integrating first over go and then 
over with the integration limits 

< < 2i?e [m - J . (37) 

For practical calculations we take Mx„.^ = on the assumption of quark-hadron duality. 



The parameters which enter in the calculation are the final quark mass m^, the pa- 
rameters a and h (or the equivalent quantities: the mean value ^ and the variance a^) of 
the distribution function and the strong coupling constant as- An important feature is 
the appearance in the decay rates of the physical 5-meson mass instead of the h quark 
mass. 

The necessity of taking into account the interplay between perturbative and nonper- 
turbative QCD effects on both 6 — > c and h ^ u spectra is discussed above. Here we 
illustrate in Fig.l that this interplay is important especially in the endpoint region of the 
b ^ u spectrum. The radiative correction calculated with the help of eq.(P^ is a smooth 
function of the electron energy up to the physical endpoint (solid curve). The other two 
curves show the quark-level perturbative correction without and with the Sudakov expo- 
nentiation [0, respectively. Radiative corrections without the Sudakov exponentiation 
run off to infinity with increasing energy as expected. The Sudakov exponentiation elimi- 
nates the singularity and gives a decay rate finite up to the quark- level endpoint = 
Beyond this value the correction is zero. For our case, the radiative correction, shown 
in Fig.l, remains finite all the way up to the physical endpoint E^, = The property 
that the averaging over a variable of the initial quark renders the radiative corrections 
finite is general and can be applied in other approaches as well. For example, radiative 
corrections to the analyses pil pSl can be easily included using our prescription. 



6 Electron Energy Spectra 



Having at our disposal a coherent treatment of perturbative and nonperturbative QCD 
effects, we study the sensitivity of the shape of the spectrum to various parameters. As 
discussed in the previous section, we use eqs.(24), (34) and (35) to compute the spectra. 
More precisely, we integrate eq.(24) over go and and then add the QCD radiative 
correction from eq.(34) to obtain the radiatively corrected electron spectrum. 



10 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



7 - 



6 r 

> 

O 5 - 




2.3 2.35 2.4 2.45 2.5 2.55 2.6 2.65 

Ee (GeV) 

Figure 1: Comparison of the radiative correction to the b ^ u electron energy spec- 
trum in the endpoint region as calculated in this paper (solid line) with the quark-level 
results of ref.[18] without (long-dashed hne) and with (short-dashed line) the Sudakov 
exponentiation for — 0.25, rrib — 4.9 GeV, — 0, a — 0.00560 and b — 0.953 . 
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0.5 1 1.5 2 2.5 3 

Ee (GeV) 

Figure 2: The shape of the electron energy spectrum from the 6 — > m inclusive semileptonic 
B meson decay in the rest frame of the B meson for various values of parameters: (1) 
m„ = 0, // = 0.93, (7^ = 0.006 (a = 0.00560, b = 0.953) (solid line); (2) m„ = 0, = 0.89, 
a2 = 0.006 (a = 0.00992,6 = 0.913) (long-dashed line); (3) m„ = 0, /x = 0.93, = 0.008 
(a = 0.00679, b — 0.960) (short-dashed line). The strong coupling constant is fixed to be 
as = 0.25 . 

The b ^ u spectrum is shown in Fig. 2 as a function of [i and u^. It is evident that the 
spectrum is much more sensitive to the mean value ^ of the distribution function than its 
variance cr^. The shape of the spectrum is insensitive to the values of the mass of the final 
quark m„ and the strong coupling constant as- The effect of the radiative corrections 
on the spectrum is shown in Fig. 3 where we have chosen two values of ag = 0.25 (solid 
curve) and cts = (dashed curve). The effect of the radiative correction is moderate. 

We study next the dependence of the shape of the electron energy spectrum for the 
b ^ c decay on various parameters. In addition to the previous parameters, the mass of 
the charm quark plays now a role. We show in Fig. 4 the electron spectrum as a function 
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Ee (GeV) 

Figure 3: The shape of the electron energy spectrum from the b ^ u inchisivc semileptonic 
B meson decay in the rest frame of the B meson for m„ = 0, /i = 0.93, — 0.006 
(a = 0.00560, b = 0.953), a, = 0.25 (sohd hue) and a, = (dashed hue) . 



13 











0.5 



1 1.5 

Ee (GeV) 



2.5 



Figure 4: The shape of the electron energy spectrum from the 6 — c inclusive semileptonic 
B meson decay in the rest frame of the B meson for various values of parameters: (1) 
m, = 1.5 GeV, /i = 0.93, = 0.006 (a = 0.00560,6 = 0.953) (solid line); (2) = 1.5 
GeV, fx = 0.89, (x^ = 0.006 (a = 0.00992,6 = 0.913) (long-dashed line); (3) m, = 1.5 
GeV, fi = 0.93, = 0.008 (a = 0.00679,6 = 0.960) (short-dashed line); (4) rric = 1.7 
GeV, = 0.93, = 0.006 (a = 0.00560,6 = 0.953) (dotted line). The strong coupling 
constant is fixed to be ct^ = 0.25 . 

of rric, /i and cr^. The shape of the spectrum is a sensitive function of the mass of the 
charm quark rric and of the mean value /i. It is rather insensitive to the variance cr^. In 
Fig. 5 we show the electron spectrum with and without radiative corrections. It is evident 
that the 6 — > c spectrum is insensitive to the value of a^. 



With the parameters determined so far we can calculate the 6 ^ c spectrum and 
compare it with the recent experimental data from the CLEO collaboration PB|. We 
present the result in Fig.6, where the theoretical curve has been boosted to the rest frame 
of the T(45') resonance. The values of a = 0.0118 and 6 = 0.931 are still consistent with 
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Figure 5: The shape of the electron energy spectrum from the b ^ c inclusive semileptonic 
B meson decay in the rest frame of the B meson for rric =1.5 GeV, = 0.93, = 0.006 
(a = 0.00560, b = 0.953), a, = 0.25 (solid line) and a, = (dashed line) . 
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Figure 6: The predicted electron energy spectrum compared with the CLEO data. The 
theoretical calculation uses = 1.61 GeV, a = 0.0118, b = 0.931, = 0.25 and 
f3 = 0.061 for the velocity of the B meson. The spectrum is normalized to the B-meson 
semileptonic branching fraction. 



the sum rules in eqs. (|3^) and (^). This is a direct calculation of the spectrum and not a 
fit. The agreement with the experimental data is good. 

The above analyses indicate the sensitivity of the spectral shape to various parameters 
and imply that a detailed fit[] to the measured spectrum can impose strong constraints 
on the mean value of the distribution function and the mass of the charm quark. This 
procedure will reduce the theoretical uncertainty in the calculation of the semileptonic 



decay width of the 5-meson [^] and improve the accuracy of the predictions for the 
h ^ u spectrum. Precise determinations of \Vcb\ and \Vub\ may be gained from inclusive 
semileptonic i?-meson decays. Finally, the same tensor structure appears in the decay 



^Such a fit should also account for detector resolution and bremsstrahlung. 



16 



B ^ J/ip + X 1^ and a universal fit of both processes would be of interest. Dedicated 
studies of the inclusive B decays will also offer more insight into the internal structure of 
hadrons. 
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